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We explore theoretically the complex dynamics and emergent behaviors of spinning spheres im¬ 
mersed in viscous fluid. The particles are coupled to each-other via the fluid in which they are 
suspended: each particle disturbs the surrounding fluid with a rot let field and that fluid flow affects 
the motion of the other particles. We notice the emergence of intricate periodic or chaotic trajecto¬ 
ries that depend on the rotors initial position and separation. The point-rotor motions confined to a 
plane bear similarities the classic 2D point-vortex dynamics. Our analyses highlight the complexity 
of the interaction between just a few rotors and suggest richer behavior in denser populations. We 
discuss how the model gives insight into more complex systems and suggest possible extensions for 
future theoretical studies. 


Collective motion of active particles has received con¬ 
siderable recent interest in the scientific community, 
spanning across many disciplines and studied using di¬ 
verse tools ranging from theoretical and computational 
to experiments. Active matter is composed of driven 
units or particles, each capable of converting stored or 
ambient energy into movement [T]. The interaction of 
such self-driven units with each-other and the surround¬ 
ing environment gives rise to intriguing spontaneous or¬ 
ganizations in their population. Biological active sys¬ 
tems include mixtures of biofilaments and motor proteins 
from cell extracts, the cytoskeleton of living cells, bacte¬ 
rial suspensions, herds of terrestrial animals such as ants, 
aquatic animals such as fish, and aerial flocks of birds . 
Nonliving active matter includes vibrated granular rods, 
an ensemble of robots, suspensions of colloids propelled 
in fluid by their surface catalytic activity when chemicals 
are added to the system, or colloids driven into motion 
by light. For more information on these examples of ac¬ 
tive systems see the recent reviews on active and driven 
matter mia and the references therein. 

In the aforementioned active systems the individual 
units are self-propellers and in most instances the gen¬ 
erated propulsion is translational. Little attention has 
been shown on active units that rotate due to internal 
or external torques, partly because such types of rotat¬ 
ing units were realized experimentally only recently. New 
experimental realizations of spinning particles have been 
reported when the particle rotation is driven by chemi¬ 
cals [3], optical tweezers or light mis], magnetic 0171 
or electric [8] fields. Increased interest in rotor systems 
generated theoretical studies exploring simple sphere ro¬ 
tors pair dynamics pun], self-assembly m, dynamics 
at interfaces [12], rheology of suspensions [niiTi]. chiral 
motor suspensions [15] or synchronization in a carpet of 
hydro dynamically coupled rotors with random intrinsic 
frequencies m- 

Other than direct particle collisions, the interactions 
in such systems are dominated by the fluid flow, i.e., at 
leading order a rotor particle is advected by the flow gen¬ 
erated by the other rotors. The fundamental interactions 
between such particles and the resulting collective dy¬ 


namics has not been investigated much, with the notable 
exception of a pair of rotors 13 uni HZ]. The pair dynamics 
has been studied on a more complex model of an active 
rotor, for example ones with an elongated shape [9], ro¬ 
tors which may be translating while also rotating m or 
Quincke rotor pairs with electro-hydrodynamical inter¬ 
actions m- The collective dynamics has been studied 
in more complex flow conditions, e.g. shear [13]. Sur¬ 
prisingly, the self-organization and coupled behavior of 
many rotor particles under no imposed external flow (i.e. 
in quiescent flow) has not been investigated. 

We study theoretically the collective dynamics of ro¬ 
tating particles and the fundamental interactions which 
lead to group behavior. For rotating particles suspended 
in viscous fluid the rot let flows dominate the interactions. 
We show how the motion of two and three particles re¬ 
sults in regular or periodic motion, as do certain sym¬ 
metrical configurations of more particles. For four or 
more rotlets, irregular or chaotic trajectories can emerge 
depending on the initial particle separations. 

Chaos is a rarity in flows dominated by viscosity - 
Stokes Equations governing the fluid flow are linear and 
obey kinematic and time reversibility - however many 
body interactions can provide a source of nonlinearity in 
the system m- For example, three Stokeslets (the sim¬ 
plest model for a sedimenting particle) can exhibit chaos 
[m |20] , whereas populations of stresslets (the simplest 
model for a micro-swimmer) can display coherent large- 
scale motions [HJ [^. Spinning particles (rotors) have 
been observed to organize in regular arrays and crystals 
[Tidiiisa] or exhibit collective directed motion [6]. The 
regular crystal structuring has been explained by the bal¬ 
ance of hydrodynamic repulsion and the magnetic attrac¬ 
tive force. 

We analyze here the collective dynamics of purely rot- 
let particles and show that unlike Stokeslets, the three 
body interaction does not lead to chaos. The similarities 
to the 2D point-vortex analogue problem are discussed. 
In the end we note the implications of pair or triplet in¬ 
teractions in the dynamics of denser rotor suspensions, 
and outline possible future research directions. 
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Model: The simplest model for a rotor is a neutrally 
buoyant sphere of radius a spinning due to the action 
of an externally imposed torque T. In the creeping flow 
limit where inertia is negligible the generated flow is 


from Eq. (§ result in planar motions only ( dzijdt = 0). 
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which describes a rotlet or couplet that decays as 1/r^ 
with distance r = |x — x^| from the rotor center (see 
m for the derivation and details). In a suspension, 
a few well-separated rotors interact with each-other di¬ 
rectly and through the fluid: each particle disturbs the 
fluid flow as it rotates and that in turn affects the other 
particles. To leading order, a rotor is simply advected by 
the flow generated by the other rotors 


dyjj 

dt 
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where Vij = |x^ — Xj| are the separation distances. 

The dynamics of a pair of rotors with aligned and same 
magnitude torques that are perpendicular to the rotors’ 
plane of motion is known um- Same-spin rotors just 
orbit around their centroid at constant speed, see Fig. 
while opposite-spin rotors translate at constant speed in a 
direction perpendicular to their separation, as illustrated 
in Fig. This form of co-operative self-propulsion 

between two rotors with different spins has been noted by 
Leoni and Liverpool [9] for non-spherical rotors and Fily 
et. al. HO] for purely rotating or rotating-and-translating 
particles. 


where tZi = Ta^/Sirfi can be positive or negative depend¬ 
ing on the rotor spin. Note that this flow field is az¬ 
imuthal, i.e. there’s no z component. The rotors remain 
in the x — y plane and the dynamics is essentially 2D. 

The (viscous) rotor dynamics Eqs. 11 looks similar 
to the the 2D (inviscid) point-vortex dynamics [25] except 
that the power of the separation in the denominator in 
Eqs. ai is 3 and not 2. 

Eqs. 11 for rotors are a Hamiltonian system 

dxi dl-L 

with the Hamiltonian 

n = 

(In 2D point-vortex dynamics it is Idy = log Hi-) 

The Hamiltonian in Eq. 1 is invariant under translation 
and rotation, and it can be derived by considering that 

tZiXi = Q = const, Kiyi = P = const, (7) 

i i 

Ki{xl + yf) = I = const. (8) 


* dt dxi ’ 


, Tij = |Xi - Xj 


(6) 



FIG. 1: (Color online) (a) A pair of same-spin rotors circu¬ 
lates around their center of mass, (b) A pair of opposite-spin 
rotors self-propels in the direction perpendicular to their sep¬ 
aration. The fluid flow is shown. Here 


Here we show that larger collections of rotating spheres 
undergo some very surprising dynamics. 

Planar motion of point rotors: Let us first analyze 
a simpler scenario where all rotors lie in the same x — y 
{z = 0) plane and all respond to a torque that aligns with 
the z-axis T = Tz. In this case the particle trajectories 


This implies that if 7^ 0 f he centroid of the point 

rotlets remains in place and can be chosen as a reference 
point in the flow. It follows from Eqs. ( [7]|^ that 

I E (9) 

^ iT^j \ i / 

is a constant of motion that is independent of the coordi¬ 
nates. This a well-known result for point-vortex motion 
found in most fluid dynamics textbooks. 

Moreover, with a straight-forward calculation, we can 
show that 



4 i^k^ijk I ^3 



( 10 ) 


where Aijk is the area of the triangle with sides 
which can be found in terms of them by 
Heron’s formula [26], “ 4j)- 

Note that the equivalent of Eq. (10) for rotors in 
the analogue 2D point-vortex dynamics involves the term 
i'^7k - ^Jk) instead. 
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Eqs. ([Tq| ) together with the two integrals of motion 
Eqs. (|^|9| specify the relative motion of N plane-bound 
rotors in terms of their separations. Moreover, the struc¬ 
ture of Eq. (10) suggests that the three rotor problem 
is an important block to understanding the dynamics for 
larger N > 3 rotors. 



FIG. 2 : (Color online) (a) An equilateral triangle configura¬ 
tion of same-spin rotors circulates around their center of mass. 
Collinear triplet of rotors: the outer ones circulate around the 
center of mass when the middle rotor has (b) same spin or (c) 
opposite spin as them. Note the different orbiting directions, 
(d) A stable equilateral triangle configuration where the ro¬ 
tor with different spin from the others circulates in the outer 
trajectory. The fluid flow is shown. In all cases = 



0 X 1 




FIG. 3: (Color online) (a) The trajectories of a triplet of 
same-spin rotors where the initial configuration (shown in 
black dots) is a horizontal perturbation of the placement of 
the middle rotor from Fig. [^b with X 3 = (1/4, 0, 0). (b) Tra¬ 
jectories of same-spin rotors when the initial configuration has 
a vertical perturbation of the placement of the middle rotors 
from Fig. with X 3 = (0,1/4, 0). (c) The time evolution 
of the set (^ 3 , r| 3 ) for three-rotor examples where the middle 
rotor initial placement is a horizontal perturbation shows as in 
(a) shows closed curves and periodicity. The points (1,1) and 
(A, 4) correspond to the stable crystal configurations in Fig. 
I^b. (d) The time evolution of the pair {rf^, for examples 
where the middle rotor initial placement is a vertical pertur¬ 
bation as in (b) also indicates periodicity. The points (1,1), 
(4,0) and (0,4) show stable crystal configurations, with the 
last two representing cases where two of the rotors overlap. 


Three rotor motion: The three-body motion is a 
classical problem in Physics and Mathematics: e.g. the 
gravitational mass analogue has been long investigated 
from Poincare himself and is one of the systems known 
to display chaos. The three-body problem in 2D vortex 
dynamics however does not display chaos [25]. Like it, 
the three rotor problem does not display chaos as it also 
belongs to a family of integrable systems. While we do 
not present a formal proof here on the integrability of the 
three rotor problem, it is possible to construct it following 
the analogue proof in 2D vortex dynamics [26]. With Eq. 
(10) for the evolution of the three separations, solving for 
the instantaneous shape of the triangle with the rotors 
as vertices now only requires elimination of variables and 
quadrature. The motion can be entirely determined if 
we have one additional equation in which one or more 
coordinates and time appear. This is not to say that 
determining the rotor trajectories is trivial. 

Two special cases of the three-rotor problem are of par¬ 
ticular interest: the identical rotors Ki = i = 1,2,3, 
or also the triplet ki = /^2 = —The initial condi¬ 
tions of the problem are very important in the resulting 


trajectories and dynamics. We investigate this question 
numerically by integrating Eqs. [Spja nd also making use 
of the constants of motion Eqs. for high accuracy. 

A stable configuration of three same-spin rotors is 
the equilateral triangle placement: if initially equispaced 
they remain so and circulate around their center of mass, 
as seen in Eig. 1^. If three rotors are initially collinear 
with one exactly at the mid-point, then the middle one 
stays in place while the outer ones circulate around it 
maintaining collinearity, as seen in Eig. [^c, no matter 
the spin of the middle rotor. The angular speed of the 
outer rotors changes however due to the influence of the 
middle rotor: circulation is faster in all cases compared 
to the two same-spin rotor case in In Eig. I^b cir¬ 
culation is counter-clockwise, in Eig. [2f it is clockwise. 

Perturbing the stable configurations of Eig. [^b by 
displacing the middle rotor horizontally (x-direction) or 
vertically (y-direction) yields an interesting triplet dy¬ 
namics, as seen in Eigs. [^b. A horizontal perturbation 
of the middle (3rd) rotor (hence initial isosceles configu¬ 
ration, yields the triplet dynamics seen in Eigj^ where 
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the trajectories are inter-mingled but nonetheless con¬ 
fined within a compact space. A vertical perturbation of 
the 3rd rotor gives an interesting dynamics as well: the 
two nearest rotors pair up and rotate around each-other 
while also circulating on the whole in a larger motion 
opposite the other rotor which traces a larger outer tra¬ 
jectory. The trajectories are periodic, and this can be 
confirmed by looking in Fig. [^d at the curves of the set 
(^ 13 :'^ 23 ) time for a variety of initial conditions such 
as those in Figs. [^b. 



FIG. 4: (Color online) (a) The regular trajectories of four 
same-spin rotors where the initial configuration (shown in 
black dots) is a rectangle, (b) Irregular (chaotic) trajecto¬ 
ries of four rotors (initial positions shown in black dots), (c) 
The time evolution of the pairs plotted in time for 

the initially-rectangular configurations as in (a). They are 
closed curves and indicate periodicity. The point (1,2) cor¬ 
responds to the stable unit square configuration, (d) The 
time evolution of the pair for the example in (b) 

suggests chaotic motion. (e,f) X 4 {t) and y 4 {t) diverge signifi¬ 
cantly from each other when the initial placement of rotor X 4 
from (b) is slightly perturbed by 10 “^ (original trajectories 
shown in black, their perturbations are shown in gray). The 
Lyapunov exponent here is positive and, since the trajectories 
of all the rotors stay in the bounded space seen in (b), this 
suggests chaotic motion. 

Four or more rotors: Just as with the two and three 
rotor cases, the four rotor problem also can have config¬ 
urations that stay stable: for example four equal-spin 
rotors that are initially placed in the vertices of a square 
will rotate around their centroid (not shown) just as the 
two or three same-spin rotor analogues. If four same-spin 
rotors are initially placed in the vertices of a rectangle, 
we see the dynamics exemplified in Fig. [^: the clos¬ 


est rotors pair-up and circulate around each-other while 
also circulating around the quartet’s centroid. The tra¬ 
jectories are periodic, and this can be ascertained by 
the method proposed by Aref and Pomphrey for four 
2D point-vortices [27]: we monitor the separations and 
record the time-dynamics of the pair {ri^^r 2 s) for initial 
ri 2 > 0 and various initial values of rl 2 = C; this results 
into a one-parameter (C) family of sections. Since the ro¬ 
tors are identical (if with same spin), any configuration 
produces 4! initial conditions with the same values of the 
integrals of motion Eqs. ([^[^. Fig. & shows examples 
of sections and was produced by perturbing in a rectan¬ 
gle configuration. The apparently smooth curves suggest 
integrability in the vicinity of the stable uniformly rotat¬ 
ing square-configuration (aspect ratio unity), and this is 
consistent with the K AM theory as well [25] . 

By contrast, a non-symmetric initial configuration of 
the four same-spin rotors yields irregular trajectories seen 
in an example in Fig. which nonetheless stay within 
a bounded space. The fourth rotor initial placement in 
this case is a perturbation of the centered equilateral tri¬ 
angle configuration. Fig. shows the Poincare section 
of the (r 43 ,r 23 ) evolution in time that results from this 
rotor initial configuration. The section shown contains 
over 10 ^ points and represents in effects a probability 
density in the phase space. This kind of random splat¬ 
ter of points, which has been produced with shorter and 
time-reversible runs as well, suggests chaotic motion, just 
like in the case of point-vortices [27[ [28] . In Fig. [^-f we 
plot the coordinates of the fourth rotor versus time in for 
the case in Fig. [^ and with a very small perturbation of 
size 10“^. The trajectories of the two cases quickly de¬ 
correlate and diverge from each-other in time, indicating 
sensibility to initial conditions. 

Effect of the third dimension: As mentioned be¬ 
fore, due to the azimuthal nature of the generated rotlet 
flow fields that dominate the interactions, the rotors re¬ 
main in the x — y plane they are initially in. This is 
also the case when the rotors are initially in different 
planes z = constant: they only move transversely since 
dzijdt = 0 due to T = Tz (see Eq. (§)• This is obviously 
a major difference from the 2D point-vortex dynamics. 

Eor an easy example, a z-direction perturbation place¬ 
ment of another rotor from the centroid of the stable 
configuration of Eig. [§ibc does not change the overall 
structure of the ensemble because of the symmetry; only 
its rotational speed is affected. The three dimensional 
dynamics of a different configuration can be more diffi¬ 
cult to predict analytically, however due to the azimuthal 
structure of the rotlet flows, the rotor trajectories remain 
in 2 : = constant planes. 

Stable rotor crystals: As we illustrated for 2-3 ro¬ 
tors in Eigs. mi for certain initial conditions and sep¬ 
arations the rotor configurations remain the same. We 
can deem them rotor crystals: rotor patterns that move 
without change of shape or size. It is possible to find crys- 
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tals for > 3 rotors by solving the equations of motions. 
For 2D point vortices many of such crystal configurations 
are well-known and well-studied; for example see the re¬ 
views [29l |30] . It is remarkable that the planar crystals 
of rotors, whose configurations can be solved from Eqs. 
®ll> are indeed similar to those of vortex crystals, even 
though the translational or rotational motion of the crys¬ 
tal is not necessarily the same. For completeness, we in¬ 
clude some planar rotor crystals as illustrations in Fig. 

Unlike the 2D vortex system however, the rotor sys¬ 
tems can have crystals that live in three dimensions. For 
example, it is easy seen that adding rotors of either spin 
at positions (0, 0, z = constant) to the crystals in Figs. 
Si would not destabilize their configurations, though 
it does affect the crystal translation or rotation. More 
elaborate 3D crystal structures can be undoubtedly be 
found analytically or numerically. 


FIG. 5: (Color online) Illustrations of some stable planar 
crystal rotor structures. Translational or rotational motion 
of the crystal is indicated with a green arrow when possible, 
as are the rotor spins (black/red for CW/CCW respectively). 






Tracer and passive particles: Just as in the (2D in- 
viscid) point-vortex case (31], it is possible to study the 
dynamics of passive particles that are inter-dispersed in 
a rotor suspension. The dynamics of tracers, or passive 
particles, though not trivial, is important for understand¬ 
ing mixing of materials at the micro-scale. Some types 
of active particles, for example micro-swimmers and bac¬ 
teria, have been explored in simulations [21] and experi¬ 
ments [35| as successful agents that can collectively mix 
a passive field at a faster rate than just by molecular dif¬ 
fusion. Rotors present an exciting possibility to achieve 
mixing. Analytically studying the dynamics of passive 
particles in a rotor ensemble is undoubtedly difficult; 
the analogue problem in 2D point-vortex dynamics, a re¬ 
stricted A/'-body problem, is complicated even for three 
vortices m- Numerical studies are however possible. 

Effect of rotor density: The dynamics of the point- 
rotors presented thus far rests on the assumptions that 
the rotors are significantly far from each-other and the 
dynamics of the flows they generate by rotation can be 
approximated by the rotlet singularity (see Eq. (§). If 
the rotor suspensions are denser (as they often are in 
experiments), then the dynamics gets far more compli¬ 


cated. The full hydrodynamics has to be accounted for, 
namely the higher order singularities in the fluid flow 
can no longer be neglected as they do modify the parti¬ 
cle translations as well as rotations at closer distances. 
Moreover direct particle collisions can occur and the lu¬ 
brication flows between the particles cannot be neglected 
either. Analytically determining the dynamics of even 
three such rotors is very difficult [24]. Numerical ap¬ 
proaches such as the Eorce Coupling Method [32] or Ac¬ 
celerated Stokesian Dynamics [33] however could be suc¬ 
cessfully employed for this task. Indeed, fluid-mediated 
interactions can significantly affect the collective rotor 
behavior [34]. Eor example, in our very recent study of 
binary mixtures of rotors with full hydrodynamical and 
lubrication interactions, we found that the motion of ro¬ 
tors is surprisingly hindered when the total rotor den¬ 
sity in increased m- Certain types of rotor clusters and 
structures, similar to those illustrated in Eig. are more 
prevalent and stable than others, but this depends on the 
suspension density [34]. The dynamics of a few rotors, as 
studied here, could prove illuminating in understanding 
more complicated structures in more complex scenarios. 

Effect of rotor shape, type and confinement: In 
experimental systems the rotor particles are not neces¬ 
sarily spherical, solid, or purely rotators. Any deviation 
from the spherical shape affects the coupled dynamics. 
In the far-held approximation non-spherical rotors ex¬ 
hibit similar dynamics to the spherical case, however a 
secondary dynamics may emerge due to the shape and 
close-range huid hows. Eor example in rotor pairs where 
each rotor consisted of two beads connected by a thin 
rod, other than the primary dynamics as in Eig. a 
weaker oscillation of each rotor is observed [9]. Note 
that in experimental systems the rotating torque does 
not necessarily remain in the same direction or constant. 
In fact, the motion of one Quincke rotor in huid by it¬ 
self is chaotic [36] if the full electro-hydrodynamics is ac¬ 
counted for. The rotational motion of a spherical colloid 
Quincke rotor gets converted into translational motion 
when they are at a surface, giving rise to the so-called 
Quincke rollers [8]. The emergence of collective motion 
in those conhned suspensions is a complex but fascinating 
phenomenon that is motivating many new studies. 

Discussion and Conclusion: We have shown that 
point rotors in viscous hows can exhibit qualitatively sim¬ 
ilar behavior to 2D point vortices, e.g., the three rotor 
motions cannot become chaotic as their dynamics is inte- 
grable and the trajectories are periodic. Our work high¬ 
lights many intriguing analogies and diherences between 
viscous rotors and inviscid vortices. The study suggests 
many potential studies stemming in such systems, e.g. 
the periodic or chaotic trajectories of rotors, the struc¬ 
ture of stable crystals, the mixing of passive particles or 
helds. Analytical and computational studies of simple 
model systems can give important insight into the dy¬ 
namics observed in more complex scenarios. 
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